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The heavy-water-moderated molten-salt reactor (HWMSR) is a newly proposed reactor concept, in which 
heavy water is adopted as the moderator and molten salt dissolved with fissile and fertile elements is used 
as the fuel. Issues arising from graphite in traditional molten-salt reactors, including the positive temperature 
coefficient and management of highly radioactive spent graphite waste, can be addressed using the HWMSR. 
Until now, research on the HWMSR has been centered on the core design and nuclear-fuel cycle to explore the 
viability of the HWMSR and its advantages in fuel utilization. However, the core safety of the HWMSR has 
not been extensively studied. Therefore, we evaluate typical accidents in a small modular HWMSR, including 
fuel-salt inlet temperature-overcooling and -overheating accidents, fuel-salt inlet flow-rate decrease, heavy- 
water inlet temperature-overcooling accidents, and heavy-water inlet mass flow-rate decrease accidents, based 
on a neutronics and thermal-hydraulics coupled code. The results demonstrated that the core maintained safety 


during the investigated accidents. 


Keywords: Heavy-water-moderated molten-salt reactor; Neutronics and thermal-hydraulics coupling; Transient analysis; Ac- 


cident analysis 


I. INTRODUCTION 


Molten-salt reactors (MSRs) use fluoride/chloride salts dis- 
solved with fissile and fertile fuel elements as fuel. During 
core operation, liquid fuel continuously circulates through the 
primary loop of the MSRs. Online reprocessing has become 
viable, through which soluble fission products (FPs) can be 
removed online and valuable heavy-metal elements can be re- 
cycled online. Online helium bubbling can also be used to re- 
move soluble fission products. Superior performance in terms 
of safety, neutronics, and fuel utilization can be achieved in 
an MSR compared with solid-fuel reactors. Most importantly, 
the intermediate isotope, Pa-233, in the reaction chain as 


Th-232 converts to U-233 (222Th 42%, 2338p, 2 CP min), 


- (27d 
233pą É OT days), 233 U) can be extracted online outside the 


core to decay to U-233. Thus, the neutron absorption of Pa- 
233 in the core can be significantly reduced. Hence, MSRs 
are considered most suitable for thorium utilization. 

The feasibility and safety of MSRs have been demonstrated 
by the successful operation of the Molten Salt Reactor Ex- 
periment (MSRE) designed and built by the Oak Ridge Na- 
tional Laboratory (ORNL), USA in the 1960s [1]. Since then, 
China (1970s), the Soviet Union (1970s), Japan (1980s), and 
the European Union (1980s) have begun to engage in technol- 
ogy development and the conceptual design of MSRs [1, 2]. 
However, most of these MSR research activities were ter- 
minated because of the lack of national-level program sup- 
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port against the background of the nuclear-industry depres- 
sion in the 1980s to 1990s. This situation improved in the 
2000s because the MSR was selected as a candidate for the six 
GEN-IV advanced reactors in 2002 [3] and many innovative 
concepts for small modular MSRs have been proposed [4]. 
Notably, in 2011, the Chinese Academy of Sciences (CAS) 
initiated the Thorium Molten Salt Reactor (TMSR) nuclear- 
energy system project with the aim of achieving the efficient 
utilization of thorium for energy production and hydrogen 
generation within the next 20- 30 years [2]. Many studies 
related to the TMSR have been conducted by the CAS [5- 
12]. 


Most of the conceptual MSR designs, as introduced above, 
belong to graphite-moderated MSRs, in which graphite with 
high temperature and corrosion resistance is used as the mod- 
erator, and fuel channels are constructed, allowing fuel salt 
to flow through the core. The results of existing research 
on graphite-moderated MSRs demonstrate their advantages 
in terms of the fuel cycle and safety; however, some stud- 
ies have shown that graphite-moderated MSRs are limited 
by the use of graphite [13]. One of these limitations is 
the positive temperature coefficient caused by graphite heat- 
ing, which leads to neutron spectrum hardening and subse- 
quently favors the fission of U-233 over the capture of Th- 
232 [14, 15]. High-performance preparation technologies are 
required for graphite to prevent the permeation of fuel and 
Xe-135 and result in a large change in the core reactivity. Fur- 
thermore, graphite must be replaced periodically because of 
neutron-irradiation damage at high operational temperatures. 
Approximately 272 tons of highly radioactive spent graphite 
must be replaced for a 1000 MWe molten-salt breeder reactor 
(MSBR) every four years, and most importantly, an effective 
approach has not been designed for the disposal of highly 
radioactive spent graphite. The use of isotopically enriched 
carbon (up to 99.9% C-12) may minimize long-term radioac- 
tive waste; however, the high cost of producing enriched car- 
bon remains an issue. Although research conducted by Li et 
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al. [16] indicates that a negative temperature-reactivity co- 
efficient may be achieved for graphite by adjusting the core 
geometric parameters, the management of highly radioactive 
graphite remains a great challenge. 


Recently, several MSR concepts have been proposed us- 
ing different moderators such as lead, ZrHy, and BeO [17] to 
address the issues associated with graphite. However, these 
alternatives have disadvantages. For example, solid modera- 
tors such as ZrH, and BeO require periodic replacement ow- 
ing to neutron irradiation. Another effective method is to re- 
place graphite with heavy water as a moderator. In 2013, re- 
searchers from the Los Alamos National Laboratory proposed 
a salt-cooled modular thorium reactor system[18]. The sys- 
tem uses heavy water as the moderator and thorium as fuel 
and is considered safer and more sustainable for thorium uti- 
lization compared to traditional solid-fuel reactors. 


Heavy water is an excellent moderator and was adopted by 
the Canada Deuterium Uranium (CANDU) reactor. It has the 
highest neutron-moderation capability (ratio of neutron mod- 
eration to absorption cross-section, 2100 [19]) among exist- 
ing moderators, approximately 12 times that of graphite. Us- 
ing heavy water as the moderator would significantly ben- 
efit both the neutron economy and nuclear fuel cycle [20]. 
Most importantly, heavy water can be efficiently purified on- 
line and recycled via a heavy-water reprocessing system, as 
in the CANDU reactor. Thus, the disposal of highly radioac- 
tive spent graphite in graphite-moderated MSRs can be ad- 
dressed. However, heavy water can extend the time of neu- 
tron moderation and diffusion compared with graphite, which 
can retard the peak of core reactivity and mitigate the fluctu- 
ation of core reactivity in an accident. Using heavy water 
as a moderator also has the additional advantage of generat- 
ing photoneutrons, which can increase the effective delayed 
neutron fraction and enhance reactor safety margins. Hence, 
safer operation can be expected for heavy-water-moderated 
MSRs (HWMSRs) [19]. 


Until now, HWMSRs have been extensively studied in 
terms of core design and fuel cycle since it was proposed 
in 2019 by Wu et al. [19]. In their study, superior perfor- 
mance of the thorium-uranium fuel cycle was observed com- 
pared with that in graphite-moderated MSRs. In 2022, Wu 
et al. conducted a study on the transition to a thorium fuel 
cycle in a 1000 MWe low- enriched uranium (LEU)-started 
HWMSR core using various scenarios, including different U- 
235 enrichment and reprocessing times (RTs). The results 
demonstrated that the transition to a thorium fuel cycle can 
be achieved by loading the LEU with 19.75 wt% U-235/U en- 
richment no lower than 5% in the HWMSR, which is impos- 
sible in a graphite-moderated MSR owing to inferior neutron 
economics [21]. The influences of the reprocessing cycle time 
and reprocessing separation efficiency on the core-actinide in- 
ventory evolution, breeding ratio (BR), and nuclear-waste ra- 
diotoxicity were also investigated in subsequent studies [21]. 
Considering the significant advantages of neutronics, the uti- 
lization of natural uranium (NU) in HWMBSR was explored by 
Wu et al. [21]. The results revealed that the fuel-utilization 
efficiency can reach 1%, which marginally surpassed that of 
the CANDU reactor. In addition, approximately 2262 kg of 
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U-233 can be produced over 20 years of operation via online 
thorium refueling and Pa-233 extraction, which is sufficient 
for starting three new thorium-fueled HWMSR cores [22]. 
In addition to large-power-scale HWMSRs, research activi- 
ties on core design and the nuclear fuel cycle in a 530-MWth 
small modular HWMSR (SM-HWMSR) were also conducted 
by Zhang et al.[23]. 

As mentioned above, the HWMSR provides considerable 
advantages in the nuclear fuel cycle, but it must be evaluated 
within the framework of safety guidelines, which have not yet 
been considered. In view of the higher realizability of small 
modular reactors (SMRs) compared with large-power-scale 
reactors, owing to the low construction cost and core-siting 
flexibility of SMRs [24], this study conducts a safety evalu- 
ation for a 530-MWth SM-HWMSR. Transients with a high 
occurrence frequency, including molten-salt overcooling and 
overheating, fuel-salt flow-rate decrease, heavy-water over- 
cooling, and heavy-water mass flow-rate decrease, are pre- 
liminarily evaluated to demonstrate the safety of HWMSRs. 


IL. OVERVIEW OF THE SM-HWMSR SYSTEM 


The SM-HWMSR analyzed in this study adopts fluoride 
molten salt dissolved with fission and fertile elements as the 
fuel and heavy water as the moderator. It has a thermal power 
of 530 MW, and the fission energy in the core is transferred by 
two loops, as shown in Fig.1[25]. The primary loop operates 
at ambient pressure and functions to safely transfer the fission 
energy to the secondary loop; it comprises the core, control- 
rod system, and moderator cooling system. Fuel-salt con- 
duits are regularly arranged in the core to separate the high- 
temperature fuel salt ( 600°C) from the low-temperature 
heavy-water moderator ( 60°C). The fission energy released 
in the fuel salt is carried out from the core by the fuel salt 
flowing through the conduits, and is then transferred to the 
secondary loop in the intermediary heat exchanger. Control 
rods installed in the core are used to regulate core reactivity 
and shut down the reactor core during accidents. Molten-salt 
drain tanks below the core allow fuel salt to be quickly dis- 
charged from the core by gravity in case of an emergency. An 
online fuel-reprocessing system was applied to remove spar- 
ingly soluble FPs, including fission gas and noble metals, and 
to recycle useful heavy-metal elements. 

To prevent the heavy water from boiling, an excellent ther- 
mal insulator, 8YSZ-50% (8 mol% Y20O3/92 mol% ZrO», sin- 
tered into porous spheres with a relative density of 50%), was 
applied to the wall of the fuel-salt conduit to minimize the 
heat transfer between the fuel salt and heavy water [26, 27]. 
According to the analysis conducted by Wu et al., a 3mm 
thickness of 8YSZ-50% is sufficient to maintain the temper- 
ature of heavy water below boiling point, which was adopted 
in this study. Moreover, two SiC layers with a thickness of 
0.5mm were placed inside and outside the 8YSZ-50% insu- 
lator layer to resist corrosion from both the fuel salt and heavy 
water [28]. A heavy-water recycling system was adopted to 
remove the deposited heat resulting from gamma and neutron 
irradiation and the heat transferred from the fuel salt. 
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Fig. 1. Outline of SM-HWMSR system. 


177 The core design of the SM-HWMSR is illustrated in 
178 Fig.2. Both the fuel salt and moderator have negative 
179 temperature-reactivity coefficients of (—1.125pcm/K) and 
180 (—3.635 pcm/K), respectively. In addition, a deep negative 
1851 void-fraction coefficient of the moderator (—60 pcm/%) is 
182 presented. The inlet temperatures for the fuel salt and heavy 
183 Water were set to 600°C and 60°C, respectively, based on 
184 the designs of the TMSR and CANDU reactors [19, 26]. A 
185 flow rate of 1.8 m? /s was used for the fuel salt with reference 
ıse to that of the MSBR [29]. The HWMSR is a liquid-fueled 
187 reactor in which fissile elements are dissolved in molten salt, 
serving as both fuel and coolant. Effective multiplication can 
be adjusted online by feeding liquid fuel, allowing it to re- 
190 Main at approximately | throughout the fuel cycle [19]. Con- 
1 sidering this, effective multiplication was not investigated in 
192 this study. Effective multiplication can be adjusted online by 
193 feeding the liquid fuel, allowing it to remain at approximately 
194 1 throughout the fuel cycle. The specific core and physical- 
195 property parameters of the fuel salt and heavy water are sum- 
198 marized in Tables | and 2, respectively. 
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Fig. 2. (a) Multi-channel model for SM-HWMSR; (b) Cross-section , 


of quarter core. 
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Table 1. Specification of main parameters for SM-HWMSR 


Parameter Value 
Core power (MWth) 530 
Core inlet temperature of fuel salt (°C) 600 
Core outlet temperature of fuel salt (°C) 650 
Inlet temperature of heavy water (°C) 60 
Outlet temperature of heavy water (°C) 62 
Inlet volume flow rate of fuel salt (kg/s) 4.7 
Heavy water velocity (m/s) 0.6 
Radius of active core (m) 1.89 
Height of active core (m) 3.78 
Number of fuel channels 267 
Thickness of SiC (mm) 0.5 
Thickness of insulator (mm) 3 

I. NEUTRONIC-THERMOHYDRAULIC COUPLING 


CALCULATIONS 


As in graphite-moderated MSRs, neutronics tightly cou- 
ples with the thermal hydraulics in SM-HWMSRs owing to 
the application of liquid fuel, in which the released fission 
energy is directly absorbed. A more rapid temperature re- 
sponse to power changes was observed in the core compared 
with solid-fuel reactors. Moreover, delayed neutron precur- 
sors (DNPs) entrained in the flowing fuel salt continuously 
circulate through the primary loop, influencing the transient 
behavior of the core during accidents. The calculation of the 
distribution of DNPs should consider the velocity of the fuel 
salt [30, 31]. The temperature distribution of heavy water is 
another issue unique to the SM-HWMSR and should be ad- 
dressed by coupling the core power and fuel-salt temperature. 
In this study, an inhouse two-dimensional dynamics code spe- 
cialized for HWMSR (HWMSR-2D) was developed based on 
the TMSR-2D, which is a dynamic calculation code designed 
for graphite-moderated MSRs [32-35], by replacing graphite 
with flowing heavy water. 


Table 2. Composition and physical property parameters for the molten salts and heavy water[19, 31] 


Parameter Value 
Fuel Salt Composition (mol%) LiF-BeF2-ThF4-UF4 
71.7% LiF - 16% BeF2 - 12.3% HNF4 
Coolant Salt Composition (second loop) (mol%) LiF-NaF-KF 
46.5%LiF - 11.5%NaF - 42% KF 
Fuel Salt Specific Heat Capacity (J/kg - K) 1985 
Fuel Salt Density (kg/m?) 3.153 - 5.92508E-4T (K) 
Fuel Salt Viscosity Coefficient (Pa - s) 8.8 
Fuel Salt Thermal Conductivity (W/m - K) 0.4461 + 0.00057’ (K) 
Boiling Point of Fuel Salt (K) 1676 
Thermal Feedback Coefficients of Fuel Salt (pcm/K) -1.125 
Heavy Water Density (kg/m?) 1.1 
Heavy Water Specific Heat Capacity (J/kg - K) 4187 
Heavy Water Thermal Conductivity (W/m - K) 0.668 
Heavy Water Viscosity Coefficient (Pa - s) 0.419 
Thermal Feedback Coefficients of Heavy Water (pcm/K) -3.635 
Allowable Temperature of Heavy Water (K) 373 
Fuel Salt Friction Coefficient 
Re < 3000 f =64/Re 
Re > 3000 f =0.3164/Re°?° 
Fuel Salt Nusselt Number 
Re < 3000 Nu = 1.62(RePrL/H)°°? 
Re > 3000 Na=(023he Pr? 
SiC Thermal Conductivity (W/m - K) 370 
Insulator Thermal Conductivity (W/m - K) 0.7514 - 0.0017T + 2E-6T” - 1E-9T® (K) 
217 A. Neutron kinetic model 233 Could be neglected because of the significantly higher veloc- 


239 ity of the neutron flux in the core. However, this should be 
Neutron kinetics is an important indicator of core safety 2° considered for DNP kinetics because the flow of the fuel salt 


oo and is defined as the behavior of spatial neutron flux over 241 Causes the DNP to flow out of the core and subsequently re- 


20 time. In this study, two-group time-dependent neutron dif- 2 duce the fraction of delayed neutrons in the core [36]. The 
21 fusion equations, in which the delayed neutrons are consid- 2° DNP balance equations in the SM-HWMSR can be expressed 


222 ered as fast-group neutrons, are employed for modeling the 74 4 follows: 
223 neutron kinetics of the SM-HWMSR: 


OC; O(UC;) 


2 
a = BD, ErP -ACi- StH 1... 
y gr 7 V Pta) + l- BYE pai + VE ya) 7 a (3) 


6 246 Where C; and A; represent the concentration and decay con- 
224 
— Medi + ` AiCi 247 stants of the 2-th group of DNPs, respectively, and U repre- 
i=1 248 sents the velocity of the fuel salt. The first term on the right- 
(1) 349 hand side of Eq. (3) represents the generation of DNPs result- 
225 250 Ing from the nuclear-fission reaction. The second and third 
1 a6 251 terms represent the decay of the DNPs and influence of the 
=e. (D2V¢2) — Si2d2 + X1-42¢1, (2) æ fuel-salt flow, respectively. 
V2 Ot 23 [n general, the neutrons in different positions of the core 
254 Make different contributions to the power, which is defined as 
25 the neutron importance. This can be calculated by solving the 
256 adjoint equations of the steady-state neutron flux [37]: 


226 


227 Where subscripts 1 and 2 represent the fast and thermal neu- 
228 tron groups, respectively, o is the neutron flux, V represents 
229 the neutron velocities, Xs and X; represent the fission and 
230 neutron absorption cross-sections, respectively, D represents 
231 the diffusion coefficient, J;_,2 indicates the transfer cross- V -(DiV¢!') —Sud! +(1- Avs PICA + gh) 
232 Section from the fast neutron group to the thermal group, v é 
233 represents the number of neutrons produced per fission, and B y ct =0 
234 B represents the total fraction of delayed neutrons. T Aa 

235 During the core operation of the SM-HWMSR, the fuel salt 

236 Continuously circulates in the primary loop at a velocity of 

27 0.8m /s. The flow effect of the fuel salt on the neutron flux 258 V - (D2V¢b) =D, + Lisod! =0 (5) 


(4) 


259 


267 


268 


274 


; 
a m 


2 
fi 2 (Esa) - Ài Ci, 


Bere = 
>s D Íy oh(r)X 


B. Thermal-hydraulics model 


Because the fuel salt flows through the fuel channels with- 
out lateral interaction in the SM-HWMSR, the parallel multi- 
channel model is adopted to evaluate the performance of the 
fuel-salt thermal hydraulics with the conservation of mass and 
momentum, expressed as follows [38]: 


ot L at 
o (Wi), 2 (WN op IW Wi 
Ot \ A; Oz \ pA?) dz 2DipA? 2pA?’ 
(9) 


where W; is the mass flow rate in the -th channel, N is the 
total number of fuel-salt channels, A; is the cross-sectional 
area of the 2-th channel, p is the density of the fuel salt, D; is 
the diameter of the 7-th channel, and Ç is the form resistance 
factor. Considering that the fuel-salt conduits are intercon- 
nected at both the inlet and outlet of the core, the pressure 
over the entire fuel-salt conduit at these two points is set to be 
the same and is considered as a boundary condition to solve 
Eqs. (8)-(9). 

Unlike the flow of fuel salt through the conduits, the 
heavy water surrounding the fuel channels connects to itself; 
therefore, lateral flow should be considered for the thermal- 
hydraulic calculation of heavy water. However, considering 
the small density change with temperature below the boiling 
point of heavy water, the lateral flow would be small; thus, 
it was not considered in this study. Lateral flow may be sig- 
nificant during severe accidents, such as fuel-channel break- 
age, which would result in the local boiling of heavy wa- 
ter. These cases will be studied in the future using a coupled 
two-phase heavy-water flow model. Considering the above, 
a single-channel thermal-hydraulics model was adopted for 
heavy-water calculations. As shown in Fig.3 the heavy wa- 
ter flows through the core with the heavy water on one side 
and molten salt on the other side of the fuel-salt channel. The 
heat transferred from the molten salt to the heavy water and 
the energy deposited in the heavy water owing to neutron and 
gamma irradiation were assumed to be removed by the heavy- 
water circulation loop. 


ee 1 ty oh, 
Ci(r) dV + (1 — — sa fy olr 


20 where x! denotes the adjoint of x. The importance of the de- 
261 layed neutrons is given by Eq. (6), which differs from that 
262 Of prompt neutrons (Eqs. (4)-(5)), because the position of the 
263 DNP changes in a timely manner, and some of the DNPs de- 
264 cay Outside the core. In addition, the effective delayed neu- 
25 tron fraction for group 2 must be weighted according to the 
26 total adjoint neutron flux, as expressed in Eq. (7): 


(r) dV 7) 
\(VX yp) gGg(r) dV 


SA 


Considering the heat transfer between the fuel salt and 
heavy water, the temperature of both can be calculated based 
on energy conservation: 


OT m(z, t) OTm(z,t)  Pm(z,t) — o(R, 2,1) 
Ot + Vin (2,1) Oz © TR pmCm TR? PmCm 
307 (10) 
OT p(z,t) OT p(z,t) _ PpalZ;t) 
308 ni (R, Z, t) 
T(R? — R?)ppep 


(11) 


where Tm is the temperature of the molten salt, Vm represents 
the velocity of the molten salt, Pm 1s the power released in the 
molten salt, Pm and cm are the density and heat capacity of 
the molten salt, respectively, R is the radius of the molten-salt 
channel, ¢(R, z, t) represents the heat flux transferred from 
the molten salt to the heavy water, T'p is the temperature of 
the heavy water, Vp is the flow velocity of the heavy water, 
Pp is the power deposited in the heavy water by fission, pp 
and cp are the density and heat capacity of the heavy water, 
respectively, and Ff, is the equivalent radius of a single lattice. 


C. Neutronics and thermal-hydraulics coupling 


The macroscopic neutron cross-sections of the material in- 
volved in the core are temperature-dependent parameters that 
should be provided before solving the neutron kinetic equa- 
tions in Section III A. In this study, they were prepared using 
the DRAGON-4 code [39], which was developed by the In- 
stitute of Nuclear Engineering of the Polytechnic School of 
Montreal. The obtained macroscopic neutron cross-sections 
were then tabulated with the temperatures of the molten 
salt and heavy water to couple the neutronics and thermal- 
hydraulics models [32]: 


~(Tim, Tp) = ABX(Tm2;, Tp2) + AV — B)S(Pin2, Tp1) 
30 + (1 — A)BE(Tmi, Tp2) + (1 — A) — B)S(Tint, Tpi) 
(12) 
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salt Vp, O 


Fig. 3. Single-channel thermal-hydraulics model for heavy water. 
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332 (13) 

As shown in Fig.4, the HWMSR-2D code comprises 
three functional modules: temperature-dependent macro- 
335 SCOpic neutron cross-section preparation, steady-state calcu- 
336 lation, and transient calculation. The temperature-dependent 
337 Macroscopic neutron cross-section module is first executed 
338 to provide the required macroscopic neutron cross-section 
339 for the steady-state calculation, in which the neutronics and 
34 thermal hydraulics are coupled by tabulated temperature- 
341 dependent macroscopic neutron cross-sections. First, a uni- 
342 form power distribution was used to obtain the initial tem- 
343 perature distributions of the fuel salt and heavy water and 
344 the flow rate in the fuel channels, which were then fed back 
345 to the neutronics calculations to update the power distribu- 
34 tion. The neutronics calculation was performed using the de- 
347 terministic method. This iterative calculation continued until 
348 the temperature of the molten salt converged, and the control- 
rod position was then iteratively adjusted to ensure that the 
core achieved criticality. Transient calculations were con- 
ducted considering the steady-state core as the initial con- 
dition. The time-dependent N-TH coupled calculation con- 
353 tinued until the preset time step was completed. The above 
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352 


354 processes were realized via programming using Fortran. The 
355 deterministic software DRAGON was used to perform cell 
neutronic calculations to provide a macroscopic cross-section 
for the HWMSR-2D. 
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D. Steady-state characteristics of SM-HWMSR 


358 


The steady state of the SM-HWMSR is the starting point 
of the transients, and should be evaluated before calculating 
the transients. Fig.5 shows the temperature distributions of 
molten salt and heavy water at the steady state with the aver- 
age flow velocity of fuel salt at 0.8m/s and that of heavy 
water at 0.6m/s. The green area in the figure represents 
the control-rod channel. The two-group neutron flux and 
power distributions are shown in Fig.6, and the flow distri- 
bution of the core is shown in Fig.7. The departure of the 
radial thermal-neutron distribution peak from the center led 
to a corresponding deviation in the core-power distribution. 
Thus, the maximum outlet temperature was located in the sec- 
ond and third groups of the fuel-salt channels. In this study, 
the transversal flow of heavy water was neglected, consid- 
ering the lower temperature of heavy water compared to its 
boiling point and the consequent marginal density difference 
across the core. In addition, the temperature difference of 
heavy water across the core was also small (< 3 °C) because 
of the excellent heat-transfer isolation of the thermal insu- 
lator. Moreover, the mass flow rate of the fuel salt within 
each fuel-salt channel varied marginally, owing to the rela- 
tively low power peak and high heat capacity of the fuel salt. 
The control rod was partially inserted into the core, which en- 
hanced the thermal neutron absorption in the core center. As 
the reactor power was more dependent on the distribution of 
thermal neutrons, the maximum power density was located in 
the second group of fuel-salt channels rather than at the center 
of the core. 
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IV. TRANSIENT ACCIDENT ANALYSIS 


387 


Transients are a type of accident that is expected to occur 
at a high frequency over the lifetime of the core and can en- 
danger the safety of the reactor; thus, they require extensive 
analysis. In the HWMSR, both the molten salt and heavy wa- 
ter have large negative temperature-reactivity coefficients, as 
presented in Table 2. A small temperature variation in either 
material may introduce a relatively large reactivity to the core 
and would cause fluctuations in core power, consequently 
impacting core safety. Additionally, malfunctioning of the 
molten-salt or heavy-water pumps would result in a reduction 
in the flow rates of the molten salt and heavy water, which 
in turn would increase the core temperature. Given the above 
considerations, this study mainly focuses on transients with a 
high occurrence, including molten-salt and heavy-water over- 
cooling accidents, a flow-rate decrease of molten salt, and 
a mass flow-rate decrease of heavy water. Considering the 
4o4 properties of structure materials and heavy water, two limita- 
aos tions, i.e., 475°C < fuel salt temperature < 700°C, and a 
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Fig. 4. Calculation flow chart of HWMSR-2D. 


4 heavy-water moderator temperature < 100°C, are considered 409 


4o7 to obtain the lower limit of transients for maintaining core 
aos Safety. 


A. Overcooling of inlet fuel salt 


40o Overcooling of the inlet fuel salt occurs as the external 
411 load suddenly increases or the flow rate in the second loop 
42 increases, which leads to the enhancement of heat transfer 


43 from the primary loop to the second loop. This transient is 
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Fig. 5. Temperature distributions of (a) molten salt and (b) heavy 
water. 


simulated by assuming that it occurs on the basis of a steady 
state and that the control rods remain still during the transient. 
Fig.8 displays the core response to a transient in which the in- 
let fuel temperature linearly decreased by 60°C during 15s, 
which refers to the study conducted by Cui et al., [32]. As il- 
lustrated in the figure, the core power and outlet temperature 
of the fuel salt gradually increased and eventually reached a 
new steady state owing to the positive reactivity introduced by 
the decrease in the fuel-salt inlet temperature. Although the 
temperature of the heavy water increased, it was only approx- 
imately 1 °C because of the excellent heat-transfer isolation 
of the thermal insulator, which was well within the safe range. 


The outlet temperature of the fuel salt decreased and then 
increased at the beginning of the transient owing to the com- 
petition between the inlet-temperature decline and the power 
increase. As the inlet temperature decreased, which intro- 
duced a positive reactivity to the core, the core power in- 
creased, leading to an increase in the outlet temperature of 
the fuel salt. However, this increase in the outlet fuel-salt 
temperature could not offset the temperature reduction result- 
ing from the decrease in the inlet fuel-salt temperature. As the 
inlet fuel-salt temperature further decreased, which indicated 
that a larger temperature increase would be imposed on the 
core power, the increase in the outlet temperature owing to the 
power increase surpassed the temperature reduction resulting 
from the inlet-temperature decrease. This resulted in a grad- 
ual increase in the outlet fuel-salt temperature. The increase 
in the fuel-salt and heavy-water temperatures, in turn, intro- 
duced a negative reactivity to the core, subsequently limiting 
the increase in core power. Therefore, the increased power 
and resulting negative reactivity owing to the temperature in- 
crease compete, leading to a new steady state with a high out- 
let fuel-salt temperature. 

Fig.9 shows the responses of the core power, fuel-salt outlet 
temperature, and heavy-water outlet temperature to transients 
in which the inlet fuel-salt temperature decreased by 20 ° C- 
100°C. As demonstrated in Fig.9(a), the case in which the 
fuel-salt inlet temperature exhibited a larger decrease resulted 
in a higher core power and longer time to achieve a steady 
state owing to the larger positive reactivity introduced to the 
reactor core. The increased power subsequently increased the 
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Fig. 6. (a), (b) Distributions of neutron flux for two-energy-group 
under steady state and (c) power-density distribution of the core. 


outlet temperature of the fuel salt, as shown in Fig.9(b). Fur- 
thermore, a decrease in the inlet fuel-salt temperature should 
not exceed 60°C to ensure core safety (outlet fuel-salt tem- 


458 perature < 700°C). The decrease in the outlet temperature 


at the beginning of the transient was owing to the abrupt 
decrease in the inlet temperature. Moreover, the increased 
temperature of the fuel salt led to an increase in the heavy- 
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Fig. 7. Mass flow-rate distribution in each group of channels. 


Orange-yellow represents fuel-salt mass flow rate, and green rep- 
resents the number of channels. 
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Fig. 8. Overcooling transient with the fuel-salt inlet temperature 
decreased by 60°C. 


42 Water temperature, but at a small scale, owing to the thermal- 
43 insulation effect of the conduit wall, as shown in Fig.9(c). 


464 B. Overheating of inlet fuel salt 


45 Inlet fuel-salt overheating occurs when the external loads 
ass suddenly decrease or the flow rate in the second loop de- 
467 creases (e.g., owing to a pipe rupture in the second loop). 

4s A 60°C increase in the fuel-salt inlet temperature triggers 
49 the frozen valve to open and discharge the fuel salt to the 
470 drain tank by gravity to ensure core safety [40]; therefore, 
471 a transient in which the inlet temperature of the fuel salt in- 
472 creases by 60°C within 15s was investigated to conserva- 
473 tively ensure core safety. Fig.10 displays the time behavior of 
474 the core power, fuel-salt outlet temperature, and heavy-water 
475 Outlet temperature during the transient. A rapid decrease in 
476 the core power was observed at the beginning owing to the 
477 rapid increase in the inlet fuel-salt temperature, which intro- 
478 duced a large negative reactivity to the core. A decrease in 
479 core power decreased the fuel-salt outlet temperature. How- 
480 ever, this decrease did not offset the increase resulting from 
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Fig. 9. (a) Power response, (b) fuel-salt outlet temperature response, 
and (c) heavy-water outlet temperature response to transients in 
which the inlet fuel salt decreased by different temperatures. 


the inlet-temperature increase at the beginning of the tran- 
sient. As the transient proceeded, the core power further de- 
creased, which caused the fuel-salt outlet temperature to grad- 
ually decrease and eventually be maintained at a value that 
was marginally higher than the inlet fuel-salt temperature. 
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Fig. 10. Overheating transient in which the fuel-salt inlet tempera- 
ture increased by 60°C. 


C. Fuel-salt flow-rate decrease 


The decrease in the flow rate of the inlet fuel salt is an an- 
ticipated transient for the SM-HWMSR and may be caused 
by a failure of or fault in the fuel-salt pump. This would in- 
crease the residence time of the fuel salt in the core, causing 
an increase in the fuel-salt temperature, and subsequently in- 
troducing negative reactivity to the core. However, it would 
mitigate the outflow of DNPs from the core, resulting in an 
increase in core power. Fig.l1 shows the responses of the 
core power, fuel-salt outlet temperature, and heavy-water out- 
let temperature to a transient in which the fuel-salt flow rate 
decreased by 60% over 15s. At the beginning of the tran- 
sient, the power increased marginally owing to an increase 
in DNPs in the core. It then decreased rapidly because of 
the significant temperature increase in the fuel salt, which 
introduced a large negative reactivity to the core. The con- 
tinuous decrease in power caused a decrease in the fuel-salt 
temperature after reaching its maximum value. Eventually, 
the core was balanced at a new steady state, with the fuel-salt 
outlet temperature reaching approximately 700°C. Periodic 
fluctuations in the DNP fraction in the core resulting from 
changes in the molten-salt flow rate led to fluctuations in both 
the core power and fuel-salt outlet temperature. Although the 
increased fuel-salt temperature increased the heat transferred 
to the heavy water, the significant decrease in core power, 
which resulted in the reduction of radiation-induced heat de- 
posits in the heavy water, led to a marginal decrease in the 
heavy-water temperature, as observed in Fig.11. 

Fig.12 shows the behaviors of power and the outlet temper- 
atures of fuel salt and heavy water during transients in which 
the inlet fuel-salt flow rate decreased by 20% to 70%. The 
fuel salt with a more decreased flow rate resulted in a longer 
residence time of the fuel salt in the core. This subsequently 
led to a higher fuel-salt temperature and introduced stronger 
negative reactivity to the core, resulting in a lower core power. 
Although a change in the DNP fraction in the core affected 
the core power, this contribution was relatively small. As a 
result, the DNP-fraction increase in the core for the transient 
with a larger decrease in fuel-salt flow rates could only lead to 
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Fig. 11. Reactor response to a transient in which the fuel-salt flow 
rate decreased by 60%. 


a greater fluctuation in power, but could not change the over- 
all decreasing trend of the core power. The decrease in the 
heavy-water temperature increased as the flow rate of the fuel 
salt decreased further because the effect of the core-power 
decrease on the heavy-water temperature was stronger than 
that of the fuel-salt temperature increase. The figures also 
show that the fuel-salt flow rate could not decrease by more 
than 60% to ensure that the core temperature was lower than 
700°C. 


D. Overcooling and overheating of inlet heavy water 


Heavy water has a high negative temperature coefficient 
of reactivity (—3.635 pcm/K). Therefore, perturbation of 
the inlet heavy-water temperature may lead to reactor acci- 
dents. Fig.13 shows a transient in which the heavy-water 
inlet temperature decreased linearly from 60°C to 40°C in 
10s. The decrease in the heavy-water inlet temperature led to 
a decrease in the bulk heavy-water temperature, which sub- 
sequently resulted in an increase in reactor power owing to 
the introduction of positive reactivity. The increased power 
then caused the molten-salt outlet temperature to increase, 
which in turn introduced negative reactivity and reduced the 
rate of the power increase. Eventually, the reactor reached a 
new steady state when the molten-salt outlet temperature ap- 
proached 700 °C. 

Fig.14 shows the behaviors of power and the outlet tem- 
peratures of fuel salt and heavy water during transients in 
which the inlet heavy-water temperature decreased by 10°C 
to 30°C. The transient with a larger decrease in the inlet 
heavy-water temperature introduced a more positive reactiv- 
ity to the core. This resulted in a higher core power, which 
subsequently led to a higher outlet temperature of the fuel 
salt. However, the increased core power and fuel-salt tem- 
perature could not offset the effect of the heavy-water tem- 
perature decrease caused by the decrease in the inlet heavy- 
water temperature. This was owing to the extremely low ther- 
mal conductivity of the insulating layer and high heat capac- 
ity of heavy water. Consequently, the outlet temperature of 
the heavy water decreased with a decrease in the inlet heavy- 
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Fig. 12. (a) Power response, (b) fuel-salt outlet temperature re- 
sponse, and (c) heavy-water outlet temperature response to transients 
in which the fuel-salt flow rate decreased by 20%-70%. 


water temperature. To maintain the core temperature below 
700 °C and ensure reactor safety, the decrease in inlet heavy- 
water temperature should not exceed 40°C. 

For transients with increasing inlet heavy-water tempera- 
ture (<100 °C at the inlet), the response of the reactor was in- 
verted to that of the transients in which the inlet heavy-water 
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Fig. 13. Reactor response to a transient in which the heavy-water 
inlet temperature decreased linearly by 20°C. 


temperature decreased and did not endanger the core safety. 


E. Heavy-water mass flow-rate decrease 


The mass flow rate of heavy water, which is similar to the 
flow rate of molten salt, plays a crucial role in reactor safety. 
If the heavy-water pump fails or a pipe ruptures in the moder- 
ator cooling system, the mass flow rate of the heavy water de- 
creases. A reduction in the mass flow rate of heavy water pro- 
longs the residence time of heavy water in the core, leading to 
an increase in the heavy-water temperature. This introduces a 
negative reactivity to the core and subsequently decreases the 
reactor power. The decreased power then results in a decrease 
in the molten-salt temperature, which introduces a positive re- 
activity to the core and mitigates the decrease in core power. 
Fig.15 shows a transient in which the heavy-water mass flow 
rate decreased linearly from 5.48 kg/s (velocity 0.6 m/s) to 
0.18kg/s in 5.8s. The heavy-water temperature increased 
significantly and subsequently introduced a negative reactiv- 
ity to the core, causing a rapid decrease in the reactor power 
and molten-salt temperature. The decreased fuel-salt temper- 
ature, in turn, triggered a positive reactivity in the core, which 
slowed the decline rate of the reactor power and molten-salt 
outlet temperature, and a new steady state was reached. Both 
the decreased core power and molten-salt temperature led to a 
marginal decrease in the heavy-water outlet temperature after 
reaching its peak of (97 °C) at approximately 130s. 

Fig.16 shows the responses of core power and the out- 
let temperatures of fuel salt and heavy water to transients in 
which the inlet heavy-water velocity decreased to 0.01 m/s — 
0.3 m/s in 5.9s-8s, respectively. The transient with a larger 
decrease in the inlet heavy water resulted in a longer residence 
time in the core. This caused a more significant increase in the 
heavy-water temperature, which resulted in larger decreases 
in the core power and fuel-salt temperature. Notably, the 
physical properties of heavy water were assumed to remain 
constant after its temperature exceeded 100°C to simplify 
the calculation. Therefore, to maintain the temperature of the 
heavy water below its boiling point and ensure the safety of 
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Fig. 14. (a) Power response, (b) fuel-salt outlet temperature re- 
sponse, and (c) heavy-water outlet temperature response to transients 
in which the heavy-water inlet temperature decreased by 10°C-— 
30°C. 


eos the core, the inlet mass flow rate of the heavy water must not 
67 decrease lower than 0.18 kg/s. 
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Fig. 15. Reactor response to a transient in which the heavy-water 
mass flow rate decreased linearly from 100% to 3.3%. 


V. CONCLUSIONS 


In 2019, Wu et al. [19] proposed a novel MSR concept 
moderated by heavy water, 1.e., the HWMSR, to address the 
limitations of conventional MSRs resulting from the appli- 
cation of graphite, such as positive temperature feedback and 
nuclear-waste management. The core design and nuclear-fuel 
cycle of HWMSR have been extensively studied; however, 
the safety of HWMSRs remains to be evaluated. Therefore, 
we conducted a preliminary safety study for an SM-HWMSR 
by analyzing typical transients specialized to HWMSRs. A 
transient-analysis code, HWMSR-2D, was developed for the 
HWMSR based on the existing code, TMSR-2D, which was 
oriented for graphite-moderated MSRs. The steady state of 
the SM-HWMSR was then calculated to provide the initial 
conditions for the transient analysis. By referring to the typ- 
ical transients of graphite-moderated MSRs and considering 
the application of heavy water, the transients induced by the 
abnormal temperatures and flow rates of fuel salt and heavy 
water, including the overcooling and overheating of inlet fuel 
salt, fuel-salt flow-rate decrease accidents, overcooling and 
overheating of inlet heavy water, and heavy-water mass flow- 
rate decrease, were chosen for analysis. In addition, we var- 
ied the transient conditions over a wide range to explore the 
safety limitations of each transient . The results indicated that 
decreases in fuel-salt temperature, fuel-salt flow rate, heavy- 
water temperature, and heavy-water mass flow rate should 
not exceed 540°C and 660°C, 60%, 40°C and 100°C, and 
3.3% for the overcooling and overheating of inlet fuel salt, 
a fuel-salt flow-rate decrease accident, the overcooling and 
overheating of inlet heavy water, and a heavy-water mass 
flow-rate decrease accident, respectively. Otherwise, severe 
accidents, including core meltdowns and heavy-water vapor 
explosions, may occur. 

However, the calculation results also demonstrated that all 
transients could reach a new stable state owing to the nega- 
tive temperature feedback, which could balance the increase 
in core power. Both the magnitude of the change and the 


s5 balance time of the new steady state were determined by the 
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extent to which the transient condition deviated from the ini- 
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Fig. 16. (a) Power response, (b) fuel-salt outlet temperature re- 
sponse, and (c) heavy-water outlet temperature response to transients 
in which the heavy-water mass flow rate decreased to 1.7%—50%. 
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656 VI. OPTIONS FOR FUTURE WORK AND 
657 IMPROVEMENTS 


es In the future, high-temperature moderators such as “LiOH 
e9 and Mg(OD)» etc can be explored to improve the concept 
660 Of MSRs. These moderators can operate at higher temper- 
e61 atures compared with heavy water, potentially enabling a 
662 more compact reactor with a reduced neutron-slowing dis- 
ses tance. Another issue is mitigating the neutron absorption of 
e64 the thermal-insulator composition of Y2O3. The use of a gas 
e65 gap filled with CO. or argon between the fuel channels and 
s Moderator would be an effective approach. 
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